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Abstract: In this paper we propose ALF gravitational instantons of types Ak and as 
models for charged particle systems. We calculate the charges of the two families. These 
are — (fc+1) for Ak, which is proposed as a model for k+1 electrons, and 2 — k for D^, which 
is proposed as a model for either a particle of charge +2 and k electrons or a proton and 
k — 1 electrons. Making use of preferred topological and metrical structures of the man- 
ifolds, namely metrically preferred representatives of middle dimension homology classes, 
we construct two different energy functionals which reproduce the Coulomb interaction 
energy for a system of charged particles. 
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1 Introduction 

In a recent paper [1], Atiyah, Manton and Schroers have proposed describing elementary 
particles in terms of 4-dimensional Riemannian manifolds having self-dual Weyl tensor. 
Electrically neutral particles are described by compact manifolds, while charged ones are 
described by non-compact manifolds which have an asymptotic circle fibration over R 3 or 
its quotient by a finite group. The class of non-compact manifolds considered includes 
some types of gravitational instanton, originally introduced by Hawking [2] as non-singular 
4-dimensional solutions of the Euclidean Einstein equations having a Riemann tensor which 
decays at infinity. Conserved quantum numbers are identified with topological properties of 
the manifold: charge is given by minus the first Chern number of the asymptotic fibration, 
baryon number by signature — but this last identification was provisional. These ideas 
were inspired by higher dimensional constructions related to the (apparently) 3-dimensional 
physics of our world, like the Kaluza-Klein circle dimension, and the derivation of approx- 
imate Skyrmions from self-dual Yang-Mills fields, i.e. 4-dimensional instantons [3]. 

In [1], concrete proposals were made for the neutron, given by CP 2 with the Fubini- 
Study metric, the proton, given by the Atiyah-Hitchin (AH) manifold [4], the electron, given 
by the Taub-NUT (TN) manifold [2]. 1 The models are purely static: time and dynamics 
are not included. The purpose of this paper is twofold: first, two infinite families of 
manifolds, namely gravitational instantons of types and D^, are presented as candidates 
for systems of charged particles. The metric tensors of these manifolds depend on a family 
of parameters which can be interpreted as the classical positions of the particles. The Af~ 
family is proposed as a model for k+\ electrons; the D^ family can describe either a particle 
of charge +2 and k electrons or a proton and h — 1 electrons. Second, two expressions for 
the total energy of the system are presented. When applied to the and D^ families, they 
reproduce the Coulomb interaction energy of the particles. Since we believe that Ak and 
Dk model systems of particles with at most one baryon, the Coulomb energy represents, 

1 In keeping with the notation used in [1], we call Atiyah-Hitchin manifold the simply connected double 
cover of the moduli space of centered Sf/(2)-monopoles of charge 2. 
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after the rest mass, the main contribution to the energy of the system. A rest mass term 
is also obtained in one of the two approaches. 

The current definition of a (non-compact) gravitational instanton, see e.g. [5], is that 
of a hyperkahler 4-manifold with curvature tensor decaying at infinity. A gravitational 
instanton is then Ricci flat and has self-dual Weyl and Riemann tensors. The definition 
is different from the original one given by Hawking, which would include solutions like 
Euclidean Schwarzschild, for which the Riemann tensor is not self-dual. Gravitational 
instantons are classified into four families [6], ALE, ALF, ALG, ALH, which have different 
volume growth at infinity [5], e.g. as r 4 in the ALE case and as r 3 in the ALF case. Since 
non-compact manifolds modelling charged particles are required to approach at infinity an 
S l fibration over IR 3 (or its quotient by a finite group), we must select the ALF case. It is 
conjectured [7] that the ALF family is exhausted by two infinite families of manifolds: of 
type Ak and of type Dk- At the topological level, the Ak, k > 1 and Dk, k > 3, manifolds 
arise as the minimal resolution of Kleinian singularities of cyclic and dihedral type. They 
retract onto a configuration of 2-spheres intersecting according to the Cartan matrix of 
the corresponding Lie algebra. Outside a compact set, the topology is that of C 2 /T where 
r = for Ak, k > 1, and T = D* k _ 2 , the binary dihedral group of order A{k — 2), for 

Dk, k > 3. Small values of k need a separate description: Aq, the TN manifold, has the 
topology of M 4 ; -Do retracts onto MP 2 and is not simply connected; D\ retracts onto S 2 
and outside a compact set has the topology of C 2 quotiented by the Z4 action generated 
by (21,-22) l— >• i(—Z2,z\); Di is the minimal resolution of sing ularities of (M 3 x S 1 )/Z 2 . 
The Ak and Dk families include the proposed models for the electron, Aq, and the proton, 
described by the AH manifold, a member of the D\ family. 

The metric of the Ak family, also known as mufti TN [8], is explicitly known: 

ds 2 = V (dr 2 + r 2 dn 2 ) + V" 1 (dij> + a) 2 (1.1) 

with 

fc+l 

jr{ lb-Pill 

where r\ > 0. Here r € [0, 00), £ [0,7r], 4> € [0, 2tt) are spherical coordinates, dfi 2 = 
d# 2 +sin 2 6 d(j) 2 and ifi is an angle having the range ip € [0, Anr]), chosen so as to avoid conical 
singularities in the metric. The point p has spherical coordinates (r,9,4>) and {pi} are the 
positions of fc+l distinct points in M 3 , the NUTs. If (x, y, z) are the Cartesian coordinates of 
the point p, related to the spherical coordinates (r, 9,<j>) in the usual way, and (xi,yi, Zi) are 
the Cartesian coordinates of pi, then ||p — Pi\\ = \/~{?c — xi) 2 + (y — yi) 2 + (z — Zj) 2 . The 
metric dr 2 + r 2 dS7 2 is the flat metric on M 3 . The one form a is locally such that da = *3dV, 
where *3 denotes the Hodge star with respect to this flat metric. Geometrically, a NUT is 
a fixed point of the action of the Killing vector d/dtp on the manifold [8]. 

The metric of the Dk family is, with the exception of Dq and D\, known only implicitly 
[9]; the D2 metric is presented as an approximation to the K3 metric in [10]. However, for 
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large r the following asymptotic approximation is known [11]: 



ds 2 = V (dr 2 + r 2 dn 2 ) + V~ l (dip + a) 2 (1.3) 

with 

where {±Pi} are 2k distinct points in M 3 , da = *3dU, 9 S [0,7r], cp G [0, 2tt), ip € [0, 2^77), 
and there is the additional Z2 identification 

(0,<£,V>) ~ (TT-e.^ + TT,-^). (1.5) 

For r — >• 00 the positions of the NUTs become irrelevant and the metric has the leading 
asymptotic form (1.3) with V given by 

V = l + ^l, (1.6) 

r 

where r = \ \p\\. Note that (1.4) is given by the superposition of a NUT of negative weight 
located at the origin and 2k NUTs at {±pi}, and is symmetric under inversion in the origin. 
Close to each NUT pi, the metric is well approximated by the TN metric. If the NUTs are 
very distant from the origin, S> 1, then eq. (1.3), with V given by eq. (1.4), reduces 
to the asymptotic form of the Dq metric. From the moduli space description of D^, see 
[7, 12], it follows that it is possible to set pk = without making the manifold singular or 
altering its topology. Doing so, eq. (1.4) becomes 

v = i-JL + yl( 1 + 1 V (1.7) 

\\p\\ f^2 \\\ P -Pi\\ \\P + Pi\\j { ' 

If the other 2{k — 1) NUTs are far from the origin, eq. (1.3), with V given by eq. (1.7), 
reduces to the asymptotic form of the AH metric. The D\ family, whose asymptotic metric 
is parametrised by pi, has been studied by Dancer [13]; for pi = one obtains the AH 
manifold. 



2 Charge 

The kinds of particle systems that can be modelled by and are strongly constrained 
by the electric charges of the manifolds, which we now calculate. All the manifolds that 
we consider have an asymptotic [/(l)-fibration over a base B. Charge Q is defined to be 
minus the first Chern number c\ of this asymptotic fibration, the integer 

Q = -oi = IJ b F, (2-1) 
where F is the 2-form field strength of a connection on the J7(l)-bundle over B. 
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For Ak, B = S 2 and the total space of the fibration is S 3 /Zfc + i if k > 1, S" 3 if A; = 0. 
For £>fc, because of the identification (1.5), B = RP 2 . If k > 3, the total space is S 3 /Dl_ 2 ; 
if A; = 2, the circle bundle over MP 2 is trivial; if k = 1, the asymptotic topology is the same 
as that of D$, but with opposite orientation. In all cases c± can be computed by considering 
a double cover D k of D k where the identification (1.5) is lifted, so that B = S 2 . Then c\ 
is half the first Chern number of the asymptotic fibration in D 

It is convenient to consider the following connection form, defined on a ?7(l)-bundle 
over S 2 , which appears in the asymptotic metrics of both A^ and D^, 

fd^ + ^Lv ont/^ 
[dtPs + {Ai) s on Efe 

where U N = S 2 \{(0, 0, -1)}, C/ s = S 2 \{(0, 0, 1)}, Vtv, V>5 G [0,2tt), = - ^- The 
gauge potentials (Ai) N , (Ai) s are given by 

(2.3) 

with (Ai)n singular for 8 = ir, (Ai)s singular for 9 = 0. The connection form ui is instead 
globally defined on the total space of the fibration. For I = 0, ujq is a flat connection on the 
trivial bundle S 2 x U(l). Since 17(1) is abelian, F[ = d(A l ) s = d(A t ) N = -{1/2) sin 9 d6 A 
d(j) is globally defined on S 2 . The first Chern number of a [/(l)-fibration over S 2 with 
connection form u>i is I, in fact 

Cl = -—[ Fi = ^- f sin#d#Ad0 = Z. (2.4) 
2vr J 5 2 4-7T J 5 2 

For both .Afc and we can identify / by rewriting the asymptotic form of the metric in 
a way that explicitly shows oji. For A^, asymptotically (1.2) becomes V ~ 1 + (k + l)rj/r, 
therefore locally a ~ (/c + l)r/ cos d0. Hence, rescaling f = r/rj, s = s/rj, ipN = tp/(2r]) + 
(k + 1)0/2 € [0, 2vr), ^5 = -0/(2??) — (fc + 1)0/2 G [0, 2vr), the metric (1.1) asymptotically 
can be written 




ds 2 ~ 



V(df 2 + f 2 dn 2 ) + 4V' 1 Uij3 N + ^±I(cos 9 - l)d<p) 2 for 9,<f>eU N 



V(df 2 + r 2 dQ 2 ) + 4V- 1 (dips + ^ (cos 9 + l)dA forfl,0€t/ s ( 2 -5) 
= V(df 2 + f W) + W' 1 (u k+1 ) 2 . 

Therefore I = k + 1 and the charge Q^ fe of Ay. is 

Q,4 fc = -ci = -(k + 1). (2.6) 

Consider now and its charge Qo k - Using (1.6), one has a ~ (/c — 2)7ycos#d0. 
This time, making the rescaling f = r/77, s = s/rj, ipN — VV 7 ? + ~~ 2)0 € [0, 27r), 
ipS = W 7 ? ~~ (k ~ 2)0 G [0, 27r), we find that the leading asymptotic form of the metric (1.3) 
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can be written 

ds 2 ~ V(df 2 + f 2 dn 2 ) + v- 1 (^ 2(fc _ 2) ) 2 , (2.7) 

with V given by (1.6). Therefore I = 2{k — 2). Dividing by 2 to pass from to Df,, we 
get the charge 

Q Dk = ---2{k-2) = 2-k. (2.8) 

The charge (2.6) of A}, is k + 1 times that of an electron (Aq) and suggests to consider 
Ak as a model for k + 1 electrons. This interpretation is supported by the form of the 
metric (1.1), which looks like Aq close to each NUT. The charge (2.8) of agrees both 
with that of a proton and k — 1 electrons and with that of a particle of charge +2 and k 
electrons. Indeed, Dk can describe both. If none of the NUTs is at the origin, then can 
be seen as the superposition of a particle of charge +2 (-Do) an d k electrons (^4o)- Because 
of the Z 2 identification of eq. (1.5), the terms in (1.4) corresponding to the Aq NUTs come 
in mirror symmetric pairs, with each pair contributing charge —1. If one pair of NUTs is 
moved to the origin, we can instead view as the superposition of a proton (AH) and 
k — 1 electrons. For both A^ and D^, Q depends only on the asymptotic topology of the 
manifold, and its value would be the same for a different configuration of NUTs giving the 
same Chern number. But this is reasonable: in a system of charged particles, the dominant 
part of the asymptotic field depends only on the total charge and not on the details of the 
configuration. 

Let us close this section with a remark on the signature r. According to the original 
proposal [1], r gives the baryon number of the manifold. For both A^ and D^, \t\ = k. The 
fact that for k > the signature of A^ does not agree with the baryon number of a system 
of k + 1 electrons suggests to drop the identification of signature with baryon number. 

3 Energy 

We now discuss how to construct energy functionals. While conserved quantum numbers of 
a particle, like charge, are associated to topological invariants of its modelling manifold, it 
seems natural to describe dynamical entities, like energy, in terms of geometrical properties 
of the manifold, like a Riemannian metric or a connection. A first possibility would be to 
consider integrals constructed out of the Riemann tensor and its contractions, but the 
fact that ALF gravitational instantons are Ricci flat and self-dual severely restricts the 
number of possible constructions, and it does not seem to be possible to get anything but 
topological quantities via this route. However, ALF gravitational instantons have non- 
trivial homology in middle dimension: H2(Mk,Z) = Z fc where stands for either A^ 
or Dk [6, 14]. Geometry allows one to select preferred representatives among the 2-cycles 
generating the homology: those of minimal area. The two energy functionals we propose 
are built from either the area or the Gaussian curvature of these preferred 2-cycles. We 
show the construction for A^_\ and D^. We work with A/ ( ._i because it has k, rather than 
k + 1, NUTs and is more easily compared to D^. It is convenient at this point to fix a 
value for the parameter r\. In order for the circles of the asymptotic fibration to have the 
same length of 2tt for both A^-i and D^, we choose r\ = 1/2 for A^-i and ij = 1 for D^. 
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Let us start with A^—i, k > 2. There is a basis of A; — 1 ordered, independent generators 
of H2{Aj i _i,'L) that are related to the simple roots of the Lie algebra ^4fc-i an d intersect 
according to its Cartan matrix. A very natural way of building representatives of these 
2-cycles as submanifolds of minimal area has been shown in [15]: consider the two NUTs 
Pi, Pi+i and orient the coordinate axes so that the line between pi and Pi+\ lies along the 
x-axis. At each point along this line there is a circle of radius 1/V(x) parametrised by 
tp. Since at the NUTs the radius collapses to zero, the surface defined by the union of all 
these circles, which we denote by S^i+i, is topologically a 2-sphere. Its (unoriented) area 
is given by 

r-2-K fXi+i 

I V~ x dil> I Vdx = 2-rr\x i+1 - Xi\ = 2-k \\p i+ i - Pi\\ , (3.1) 

a constant factor times the Euclidean distance between the two NUTs. If instead of a line 
we had taken any other curve connecting pi and Pi+i, the area would have been the same 
constant factor times the Euclidean length of the curve, hence the 2-cycle built above is of 
minimal area. This construction can evidently be done for any pair of distinct NUTs pi, 
Pj, giving the minimal 2-cycle Sij. While {S± : 2 , • • • , Sk-i : k} is a basis for #2(^-1, ah 
pairs of distinct NUTs play an equal role in A^_i and it would be unnatural to consider 
only the basis 2-cycles. The natural choice is to consider instead all the minimal 2-cycles 
connecting pairs of distinct NUTs. Such a choice is invariant under the Weyl group of the 
Lie algebra Ak-i, the symmetric group Sk- The sets {pi}, i = 1, . . . , k, of NUT positions 
and {e^}, i = 1, . . . , k, of canonical basis vectors of M fc have the same cardinality and so a 
one-to-one correspondence pi can be established between them. A possible choice for 
the set of all roots of A^-i is {±(ej — ej), i 7^ j = 1, . . . , k}. The root — ej is associated 
with Sij and its negative with the same 2-cycle taken with the opposite orientation, Sjj. 
Our first proposal for an energy functional is simply to take tt times the sum of the inverse 
areas of the minimal 2-cycles S%j: 

< } = T i1 ~ TT' (3-2) 

Ak - x 4^ 1 \Pi -Pi\\ 



Our second proposal for an energy functional still involves the 2-cycles Si j. The 
endpoints of these 2-cycles, the NUTs pi and pj, are geometrically preferred points. It is 
then interesting to calculate the Gaussian curvature Kij of Sij at the points pi and pj. 
Without loss of generality, we can relabel these points so that i = 1, j = 2 and orient 
the axes so that p\ = (xi,0, 0), P2 = (x2,0,0). The metric on 5^2 is that of a surface of 
revolution, ds 2 = Vdx 2 + V~ 1 dtp 2 . Its Gaussian curvature Ki^ is independent of ip and is 
given by, see e.g. [16], 

Ki,2 - — 2 (3.3) 

In calculating (3.3) at p\ we take advantage of the fact that many terms are zero. Introduce 
the shorthand notation di = \\p — Pi\\, do = \\p\\, dij = \\pi — Pj\\ and group the terms in 
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V 1 under a common denominator: 



V 



f 



f + g' 



f = U d 



1 k 

-*Y^d 1 ...d l ...d k 



1=1 



i=i 



where terms with a hat are to be omitted. Then 
1 1 



2(f + g) 3 

In the limit p — > p\ one has 



(f xx g - fgxx)(f + g) - 2(/i + g x ){fxg - fg x )\ ■ 



d, 



dn for I ^ 1 



I otherwise 
Therefore in this limit 



dch 
dx 



for/^1 
1 otherwise 



dx 2 



yf+± 



(3.4) 



(3.5) 



for I ^ 1 

otherwise 
(3.6) 



/ ->• 0, / a ->■ d 2 l • • • dfcl, /ra -> 2(d 2 . . . dfc)* 



pi 



1 1 
5 -> -^21 • • • d k x, g x -> -(d 2 ...d fc ) a 



The expression for Ifi 2 simplifies to 
1 



Pi j=2 



■ ■ ■ dn . . . dki- 



2<7 
1 



\ fxxg ^f X \fx gx. 



2{d 21 ...d kl y 



2(d 2 i ■ ■ ■ d k i) x ■ -dix . . . (4i - 2(d 2 i . . . d kl ) 2 



hi ( g^ 21 • • • d kl)x + ^ ( d 31 • • • 4l + ^21^41 • • • d k i + 



1 1 1 



2d 21 2d 3 i 
More generally, for any 2-cycle Sij, 



2d k i 



»J Pi) I 2d u 2d 2 % 2du 



2d kl 



(3.7) 



(3.8) 



(3.9) 



Note that (3.9) does not involve the point pj in any particular way, so that Kij(pi) = 
K i: i(pi) for any I ^ i. Therefore we define K{pi) = K^iipi) for any 2-cycle S^i with 
/ ^ i. Summing over all the NUTs, dividing by 4, and reverting to the original notation 
dij = \\Pi- Pj 1 1 we get 



E 



^ = l y K{ ) = k+ y _L_ 



i=l 



i<i= 



1 \\pi-pj\l 



(3.10) 
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Expressions (3.2) and (3.10) are very similar. In both cases the Coulomb interaction 
energy of a system of k particles having the same charge is reproduced and the construction, 
which involves only pairs of electrons, does not run into self-energy problems. The only 
difference between 1 an d Ej^ : is that the latter contains an additive constant equal 
to the number of electrons. It seems natural to relate the length of the asymptotic circles 
to the classical electron radius r e = e 2 /(m e c 2 ), where e and m e are the charge and rest 
mass of the electron. Equation (3.10) is dimensionless. In dimensional form, with rj = r e /2, 
and multiplying by an overall constant factor m e c 2 , we get 



2 

! \\Pi~Pj\ 



E%l_ 1 =km e c 2 + £ n^^TiT, (3.H) 



the sum of the rest masses and the interaction energy of k electrons. 

Consider now the family, k > 2. Since our analysis is based on eq. (1.3), which 
holds only asymptotically, we expect the energy functionals that we derive to be accurate 
if all the pi are far from the origin. In order to construct Ep, we need to sum the inverse 
areas of the 2-cycles corresponding to all the roots of the Dk Lie algebra. A possible choice 
of roots is {±(ej ± e^), i < j = 1, . . . , k}. 2 As before, we have the correspondence pi -H- e$, 
hence the 2-cycles S±i j ±j, connecting the NUTs to the NUTs ±py, correspond to the 
roots ±ej ± ej. Multiplying by ir and summing the inverse areas, we obtain 

Dk ^VllPi-Pill \\Pi+Pj\\J 

As mentioned in section 2, we can relate Dk to either k — 1 electrons and a proton or 
k electrons and a particle of charge +2. However, eq. (3.12) contains only interaction 
terms of electron-electron type, i.e. involving the distances ||pj±Pj-||, i ^ j, and no 
interactions of proton-electron type, which involve the distances \ \pi — 0|| and should come 
with the opposite sign. While it is possible to construct additional 2-cycles corresponding 
to the proton-electron interactions, it is difficult to get the right numerical factors. Further 
investigation of the area of these 2-cycles would be worthwhile, but for the moment we opt 
not to consider E^ as a suitable energy functional for D^. 

Let us, instead, construct E K D ' . As before, the Gaussian curvature of a 2-cycle Sij can 
be calculated using eq. (3.3). We set Pi+k = ~Pu i = 1, ■ ■ ■ ,k, and 

y- 1 = — — / = TTdj, g = --r, h = - V^dodi . . .di . . .d 2k . (3.13) 

f + g + h ^ do 4^ 

Proceeding similarly as we did for Ak-i, we find that the Gaussian curvature of the 2-cycle 



2 The cases fc = 2, 3 are degenerate. For k = 3, the Lie algebra D3 is isomorphic to A3. The equivalence 
can be checked by verifying that the Cartan matrix of D3 associated to the ordered set of simple roots 
{e2 — 63, ei — e2, e2 + 63} is equal to the Cartan matrix of A3. Similarly one can show that D2 is isomorphic 
to A± x Ai. 
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i 7^ j, i^j + k, j^i + k, in the limit p — >• pi is 

^(^=4 ( l_-l T+ l J2 ( 1 + 1 H. (3.14) 

In order to get (3.14) we have used dij = \\pi — Pj\\ if either j,i < k or j,i > k, = 
lb» otherwise, and doj = \\Pj\\- Eq. (3.14) does not involve in any particular way, 

so we can again define K{p{) = Kij{pi) for any 2-cycle Sij with i^j,i^j + k, j^i + k. 
Note that K(pi) = K(pi + k). Since pi is identified with pj+fc, we sum K(pi)/A over pi for 
i = 1, . . . , k only, obtaining 

< 4 1 km = * - 1 ^ + ,£ + ra) • (3 - 15) 

If — <C lbi||; Ibill; the term 1/ +jPj|| is negligible, and (3.15) reduces to an 
additive constant, equal to the number of electrons, plus the Coulomb interaction energy 
of k electrons and a particle of charge +2. If all the pi are far from the origin, corrections 

(2) 

to Ejj due to the different behaviour of the exact metric near the origin are small and we 
expect them to be related to the rest mass of the positively charged particle. In dimensional 
form, with rj = r e , and multiplying by an overall constant factor m e c 2 , 

«S=^-gS^(Fi^ira)' (3 ' 16) 



Let us now see what happens to E D ^ in the limit p^ — > 0. We cannot take this 
limit directly in (3.15), but we can calculate the Gaussian curvature of the 2-cycles S±i t ±j 
using the metric (1.3) with V as in (1.7). Nothing particular happens to the 2-cycles 
{S±i t ±j,i,j = 1, . . . , k — 1, i 7^ j}. The 2-cycles {5±t,±ifci i = I, ■ ■ ■ , k — 1} now connect the 
points {±Pi} to the origin. Their geometry near the origin is not described accurately by 
the metric (1.3), but the description is accurate near pi. The Gaussian curvature at the 
point pi / of any 2-cycle Sij, i^j,i^j + k, j^i + k, calculated using eq. (3.3), with 
V as in (1.7), is 

KM ^{'-m + 5 jt (ra + isrbsi)) ■ (3 - 17) 

Therefore, neglecting the contribution of the origin, we have 
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In dimensional form, 



E® = (k — l)m e c 2 -Y-^-+Y ( e * + e ' ) . (3.19) 

We expect the contribution of the origin to modify the rest mass term; the interaction 
terms in (3.19) can be recognised as the Coulomb interaction energy of a proton and k — 1 
electrons. 

4 Conclusions 

In [1], self-dual 4-dimensional manifolds were used in an attempt to provide a geometri- 
cal description of matter. Their various topologies are potentially related to the particle 
quantum numbers, like baryon and lepton number, and electric charge. Their geometry is 
potentially related to particle energies and interactions. Many issues were left open, for 
example dynamics was not considered and the number of concrete examples proposed was 
relatively small. 

In this paper we have started filling some of the gaps. The class of concrete examples 
has been extended by two infinite families, the and ALF gravitational instantons. 
The former has been proposed as a model for fc+l electrons, while the latter for a positively 
charged particle and some electrons. The positively charged particle can be a proton, in 
which case there are k — 1 electrons, or a particle of charge +2, in which case there are k 
electrons. As a first step towards dynamics we have constructed energy functionals which 
reproduce the appropriate Coulomb interaction energy for the particle systems considered. 
Our formulae are precise for the Ay. manifolds, whose metric is known exactly, but only 
approximate for the manifolds, whose metric is only accurately known away from the 
origin. 

Our constructions have hinted at the important role of 2-dimensional substructures in 
this geometrical approach. Finally, this work supports the identification of charge with the 
first Chern number of the asymptotic fibration but suggests a rethinking of the identification 
of baryon number with signature. 
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